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Most recently, both BaBar and Belle experiments found evidences of neutral D mixing. In this 
paper, we discuss the constraints on the strong phase difference in D'^ Ktt decay from the 
measurements of the mixing parameters, y' , ycp and x at the B factories. With CP tag technique at 
i/'(3770) peak, the extraction of the strong phase difference at BES-III are discussed. The sensitivity 
of the measurement of the mixing parameter y is estimated in BES-III experiment at ■!/)(3770) peak. 
Finally, we also make an estimate on the measurements of the mixing rate Rm. 
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Due to the smallness of AC — amplitude in the 
Standard Model (SM), D'^ — mixing offers a unique 
opportunity to probe flavor-changing interactions which 
may be generated by new physics. The recent measure- 
ments from BaBar and Belle experiments indicate that 
the D'^ — mixing may exist 0, At the B fac- 
tories, the decay time information can be used to ex- 
tract the neutral D mixing parameters. At t = the 
only term in the amplitude is the direct doubly-Cabibbo- 
suppressed (DCS) mode D° isT+Tr", but for i > 
mixing may contribute through the sequence 
Z)0 K^TT" , where the second stage is Cabibbo 
favored (CF). The interference of this term with the DCS 
contribution involves the lifetime and mass differences 
of the neutral D mass eigenstates, as well as the final- 
state strong phase difference 6kt!- between the CF and 
the DCS decay amplitudes. This interference plays a 
key role in the measurement of the mixing parameters at 
time-dependent measurements. 

With the assumption of CPT invariance, the mass 
eigenstates of 13° — 1)° system are \Di) = p\D'-*) + q\D'^) 

and \D2) = p\D'^)—q\D'^) with eigenvalues — toi — -Fi 

and fi2 = rn2 — 2^^' ''^^P^'^^i^^^yi where the mi and Fi 
(m2 and F2) are the mass and width of Di {D2). For 
the method of detecting — mixing involving the 
— > Kt: decay mentioned above, in order to separate 
the DCS decay from the mixing signal, one must study 
the time-dependent decay rate. The proper-time evolu- 
tion of the particle states l-Cpjjyg(i)) and l-C'pjjyg(t)) are 
given by 



with definitions 
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_ ■mi+m2 . 
m = , /\m — m2 — mi, 

F = ^^i^, AF-F2-F1, 
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Note the sign of Am and AF is to be determined by 
experiments. 

In practice, one define the following mixing parameters 



_ Am _ AF 

^ = ~' ^= 2r- 
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The time-dependent decay amplitudes for ^p}jyg(0 " 
if+TT^ and ^p}jyg(i) K~tt'^ are described as 



^^K+.-[Ag+(i)-5-W], (5) 
P 



{K-7T+\n\Dly^y^{t)) - g^{t)AK-.+ -^9-it)AK-.+ 



= ^A^-^+[Xg+it)-g-it)], (6) 

where A^^- = {K+tt-\H\D") , Ak+^- 
{K+TT-\n\D^), Ak-^+_ = {K-TT+\n\D^), and 
Ak-.^+ = {K-T:+\n\D"). Here, A and A are de- 
fined as: 
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From Eqs. ([5|) and jS]), one can derive the general ex- 
(2) pression for the time-dependent decay rate, in agreement 
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with 1,3: 

dtAf 



[(|A|2 + l)cosh(j/n) + 
(|A|2 - l)cos{xVt) + 
27^e(A)sinh(y^^) + 
22:m(A)sin(a;rt)] (9) 



dtN 



\A 



K--IT+ I 



[(|Ap + l)cosh(yrt) + 
(|Ap - l)cos(xR) + 
27^e(A)sinh(^/^^) + 
2Jm(A)sin(a;rt)] (10) 

where A/" is a common normahzation factor. In order to 
simplify the above formula, we make the following defi- 
nition: 



We can characterize the CP violation in the mixing 
amplitude, the decay amplitude, and the interference 
between amplitudes with and without mixing, by real- 
valued parameters Am-, Ao, and <j) as in Ref d i]. 
In the limit of CP conservation. Am, Ad and are 
all zero. Am — means no CP violation in mixing, 
namely, \qlp\ = 1; Ab = means no CP violation in 
decay, for this case, r = r' = Bd = \Ak-'k+ / Ak'-k+\^ — 
{Ak+tt- / Ak+tt- = means no CP violation in the 
interference between decay and mixing. 

In experimental searches, one can define CF decay as 
right-sign (RS) and DCS decay or via mixing followed 
by a CF decay as wrong-sign (WS). Here, we define the 
ratio of WS to RS decays as for Z?°: 



Kit) 



and for Z?°: 



Kit) 



dWxe-'^W X 2\Ak~^+\^'' 
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^ = {l + AM)e-'^, 
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where (3 is the weak phase in mixing and Am is a real- 
valued parameter which indicates the magnitude of CP 
violation in the mixing. For / = K^n'^ final state, we 
define 
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where r' and a' (r and a) are the ratio and relative phase 
of the DCS decay rate and the CF decay rate. Then, A 
and A can be parameterized as 



A = -V? 



1 



1 + Am 



A = -V^(l + AM)e-*("+''). 



(13) 
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In order to demonstrate the CP violation in decay, we 

+ Ad) and ^f? = ^/r^—^—. 

t + Ad 



(15) 



Thus, Eqs. (fTSj) and (fT4|) can be expressed as 



1 + Am 



1 + Am ,(s+4>) 
""i + Ad 



(16) 



where S = — - — is the averaged phase difference be- 



tween DCS and CF processes, and 



■p. 



Taking into account that |A|, |A| <C 1 and x, y <^ 1, 
keeping terms up to order x^, and Rd in the ex- 
pressions, neglecting CP violation in mixing, decay and 
the interference between decay with and without mixing 
(Am = 0, Ad = 0, and = 0), expanding the time- 
dependent for xt, yt < r~^, combing Eqs. ([9]) and pO]). 
we can write Eqs. ([T7| and (fT8|) as 



R{t) = R{t) = Rd + ^/R^y'Tt 



.(rt)2, (19) 



where 



x' — xcosS + ysinS, 
y' — —xsinS + ycosS. 



(20) 
(21) 



In the limit of SU(3) symmetry, Aj^+t^~ and Aj^+t^- 
{Ax-T^+ and A^-^+) are simply related by CKM fac- 
tors, {VcdV*jVcsV*^)Aii+^- 0. In particular, 
Ak+tt- ^-iid Ax+TT- have the same strong phase, leading 
to a' = a = in Eq. ((T2]). But the SU(3) symmetry 
is broken according to the recent precise measurements 
from the B factories, the ratio Q: 
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(22) 



is unity in the SU(3) symmetry limit. But, the world 
average for this ratio is 



Uerp = 1.21 ±0.03, 



(23) 



computed from the individual measurements using the 
standard method of Ref. Since the SU(3) is bro- 

ken in 13 Ktt decays at the level of 20%, in which 
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case the strong phase 6 should be non-zero. Recently, a 
time-dependent analysis in D Kit has been performed 
based on 384 fb~^ luminosity at T(45') By assuming 
CP conservation, they obtained the following neutral D 
mixing results 



Rd — 



(3.03 ±0.16 ±0.10) X 10"^ 
(-0.22 ±0.30 ±0.21) X 10"^ 
(9.7 ±4.4 ±3.1) X 10"^ 



(24) 



on the DP are model-dependent, the results suffer from 
large uncertainties from the DP model. In this analysis, 
they see a significance of 2.4 standard deviations from 
non-mixing. Here, we will use the value of x measured in 
the DP analysis for further discussion. As shown in Eq. 
(|2ip . once y, y' and x are known, it is straightforward to 
extract the strong phase difference between DCS and CF 
decay in D° Ktt decay. If taking the measured central 
values of x, ycp{~ y) , and y' as input parameters, we 
found two-fold solutions for tanJ as below: 



tan(5 = 0.35 ± 0.63, or -7.14 ±29.13, 



TABLE I: Experimental results used in the paper. Only one 
error is quoted, we have combined in quadrature statistical 
and systematic contributions. 



Parameter 


BaBar (xlO"'') 


Belle(xl0"-') 


Technique 


X 

y' 

Rd 

ycp 

X 

y 


-0.22 ±0.37 [IJ 

9.7 ±5.4 [IJ 
3.03 ±0.19 [Ij 


0.18l»-JJ [8] 

o.6it.^ m 

3.64 ±0.17 [8j 
13.1 ±4.1 [2] 
8.0 ± 3.4 [9] 
3.3 ± 2.8 [9] 


Ktt 
Ktt 
Ktt 

K + K-, TT + TT- 

KsTT^TT^ 

Kstt'^tt^ 



The result is inconsistent with the no-mixing hypoth- 
esis with a significance of 3.9 standard deviations. The 
results from BaBar and Belle are in agreement within 2 
standard deviation on the exact analysis of y' measure- 
ment by using D — > Ktt as listed in Table [J As indicated 
in Eq. ([23| , the strong phase 5 should be non-zero due to 
the SU(3) violation. One has to know the strong phase 
difference exactly in order to extract the direct mixing 
parameters, x and y as defined in Eqs. (g]). However, 
at the B factory, it is hard to do that with a modcl- 
indepcndcnt way 0, [13] ■ In order to extract the strong 
phase 5 we need data near the DD threshold to do a CP 
tag as discussed in Ref. [3] . Here, we would like to figure 
out the possible physics solution of the strong phase 5 by 
using the recent results from the B factories with differ- 
ent decay modes, so that we can have an idea about the 
sensitivity to measure the strong phase at the BES-III 
project. 

In Ref 0, Belle collaboration also reported the result 

_ T{D°^K+r 



of ycp = 



1, where fcp = K^K and 



ycp = (13.1 ±3.2 ±2.5) x 10" 



(25) 



The result is about 3.2(t significant deviation from zero 
(non-mixing). In the limit of CP symmetry, ycp = U tlHj 



[l3 |. In the decay of — > Kstt^tt , Belle experiment 



has done a Dalitz plot (DP) analysis [9j, 
the direct mixing parameters x and y as 



they obtained 



X = (8.0 ± 3.4) X 10-^ y = (3.3 ± 2.8) x 10~^ (26) 

where the error includes both statistic and systematic un- 
certainties. Since the parameterizations of the resonances 



(27) 



which are corresponding to (19 ± 32)° and (-82° ± 30)°, 
respectively. 

At ■0(3770) peak, to extract the mixing parameter y, 
one can make use of rates for exclusive £>°£)° combina- 
tion, where both the Z?° final states are specified (known 
as double tags or DT), as well as inclusive rates, where 
either the 7?° or £)° is identified and the other D° de- 
cays generically (known as single tags or ST) [13]. With 



the DT tag technique [ij, ll5|, one can fully consider 
the quantum correlation in C = — 1 and C = +1 Z)°£'° 
pairs produced in the reaction e^e~ D°£'°(n7r°) and 
e+e" D°W-i{mT°) [ll,[ll,[l3, respectively. 



For the ST, in the limit of CP conservation, the rate 
of _D° decays into a CP eigenstate is given as 131 ]: 



F/,, ^ T{D^ ^ /,) = 2Al [1 - in] , 



(28) 



where is a CP eigenstate with eigenvalue 7y = ±1, and 
Af^ = \{fjf\H\D^)\ is the real-valued decay ampfitude. 

For the DT case, Gronau et. al. [7| and Xing ^ 
have considered time-integrated decays into correlated 
pairs of states, including the effects of non-zero final state 
phase difference. As discussed in Ref. [7|l, the rate of 
(£,o^O)C=-i ^ {l±X){fn) is described as Q: 



(29) 



where Ai±x = ](^*-'^]^]^°) ] is real-valued amplitude for 
semileptonic decays, here, we neglect y^ term since y <C 
1. _ 

For C = — 1 initial D'^D'^ state, y can be expressed in 
term of the ratios of DT rates and the double ratios of 
ST rates to DT rates El: 



1 f^iJ+^f- ^i-J-'^U 



(30) 



For a small y, its error, A(?;), is approximately 1/ \/Ni±x, 
where Ni±x is the total number of {l^X) events tagged 
with CP-even and CP-odd eigenstates. The num- 
ber Ni±x of CP tagged events is related to the to- 
tal number of D°_D° pairs N{D"D") through Ni±x ~ 
NiD"D")[Bn{D° ^l±+X)x Bn{D" -> /±) x etag] ~ 
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1.5 X 10~'^iV(D°D°), here we take the branching ratio- 
times-efficiency factor {BTZ{D° f±) x etag) for tagging 
CP eigenstates is about 1.1% (the total branching ratio 
into CP eigenstates is larger than about 5% Q)- We find 



A(y) 



±26 



N{DODO) 



±0.003. 



(31) 



If we take the central value of y from the measurement 
of ycp at Belle experiment 2], thus, at BES-111 exper- 
iment pJl], with 20/6~^ data at iplSTTO) peak, the sig- 
nificance of the measurement of y could be around 4.3 ct 
deviation from zero. 

We can also take advantage of the coherence of the 
mesons produced at the ijj{3770) peak to extract the 
strong phase difference S between DCS and CF decay am- 
plitudes that appears in the time-dependent mixing mea- 
surement in Eq. pH)) 0, [l3| . Because the CP properties 
of the final states produced in the decay of the ^^(3770) 
are anti-correlated [H, |T3| , one 1?° state decaying into a 
final state with definite CP properties immediately iden- 
tifies or tags the CP properties of the other side. As 
discussed in Ref. 0, the process of one D° decaying to 
K^TT^, while the other Z)° decaying to a CP eigenstate 
frj can be described as 



(32) 



where A = \{K-7r+\H\D")\ and A/„ = |(/^|7^|i:>") | are 
the real-valued decay amplitudes, and we have neglected 
the terms in Eq. ([5^ . In order to estimate the total 
sample of events needed to perform a useful measurement 
of S, one defined 0, [13] an asymmetry 



(33) 



where Vktv-j^ is defined in Eq. ([5^ . which is the rates for 
the V'(3770) D'^D^ configuration to decay into fiavor 
eigenstates and a CP-eigenstates f±. Eq. ([5^ implies a 
small asymmetry, A = 2\/R]jcosS. For a small asymme- 
try, a general result is that its error is approximately 
l/V-^K-TT+i where Nk-tt+ is the total number of events 
tagged with CP-even and CP-odd eigenstates. Thus one 
obtained 



A(cos5) 



1 



(34) 



The expected number Nx-tt+ of CP-tagged events 
can be connected to the total number of D^D^ pairs 
N{D°D°) through Nk-^+ « N{D'^D°)Bn{D° 
K-TT+) X Bn{D° ^ /±) X etag ~ 4.2 X 10-^N{D°D°) 0, 
here, as in Ref 7i|, we take the branching ratio-times- 
efficiency factor 67^(P)0 ^ /±) x etag = 1-1%. With 




FIG. 1: Illustrative plot of the expected error (AS) of the 
strong phase with various central values of cos5. The expected 
error of cosS is 0.04 by ssuming 20/6~^ data at V'(3770) peak 
at BES-III. The two asterisks correspond to S — 19" and 
—82°, respectively. 



the measured Rd = (3.03 ± 0.19) x 10"^ and BTZ{D" 
K-t:+) = 3.8% 4], one found 0] 



A(cos(5) 



±444 



N{D°D^) 



(35) 



At BESlll, about 72 x 10^ D°D^ pairs can be collected 
with 4 years' running. If considering both K^tt^ and 
K^t:^ final states, we thus estimate that one may be 
able to reach an accuracy of about 0.04 for cosS. Fig- 
ure [1] shows the expected error of the strong phase 5 
with various central values of cos5. With the expected 
A(cos(5) = ±0.04, the sensitivity of the strong phase 
varies with the physical value of cosS. For (5 = 19° and 
—82°, the expected error could be A((5) = ±8.7° and 
±2.9°, respectively. 

By combing the measurements of x in Kstttt and 

ycp from Belle, one can obtain Rm = (1.18±0.6) x 10"''. 
At the -0(3770) peak, D^D^ pair are produced in a state 
that is quantum- mechanically coherent 3, l3| • This en- 
ables simple new method to measure -D° mixing param- 
eters in a way similar proposed in Ref. At BES-III, 
the measurement of Rm can be per formed unambigu- 
ously with the following reactions [16|: 



{i) e^e^ 


V(3770) - 






(a) 


V^(3770) - 






[Hi) 


D-D*+ - 


{K+n- 





(36) 
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Reaction (i) in Eq. ((36)) can be normalized to D^D^ 
(iir~7r+)(i4r+7r~), the following time-integrated ratio is 
obtained by neglecting CP violation: 



N[{K-tt+){K-t:+)] ^x^ + 
N[(K-t:+)[K+t:-)] ^ 2 



= Rm. 



(37) 



For the case of semileptonic decay, as {ii) in Eq. ([M)) . we 
have 



Rm, 



(38) 



The observation of reaction (j) would be definite evi- 
dence for the existence of — D° mixing since the final 
state {K^Tr^){K^Tr^) can not be produced from DCS 
decay due to quantum statistics [la . [l7j . In particular, 
the initial D^D^ pair is in an odd eigenstate of C which 
will preclude, in the absence of mixing between the 
and over time, the formation of the symmetric state 
required by Bose statistics if the decays are to be the 
same final state. This final state is also very appealing 
experimentally, because it involves a two-body decay of 
both charm mesons, with energetic charged particles in 
the final state that form an overconstrained system. Par- 
ticle identification is crucial in this measurement because 
if both the kaon and pion are misidentified in one of the 
two ZJ-meson decays in the event, it becomes impossi- 
ble to discern whether mixing has occurred. At BESIII, 
where the data sample is expected to be 20 fb~^ inte- 
grated luminosity at V'(3770) peak, the limit will be 10""* 
at 95% C.L. for Rm, but only if the particle identification 
capabilities are adequate. 

Reactions {ii) and {Hi) offer unambiguous evidence 
for the mixing because the mixing is searched for in the 
semileptonic decays for which there are no DCS decays. 
Of course since the time-evolution is not measured, obser- 
vation of Reactions {ii) and {Hi) actually would indicate 
the violation of the selection rule relating the change in 
charm to the change in lep tonic charge which holds true 
in the standard model [16| . 

In Table [ill the sensitivity for Rm measurements in 
different decay modes are estimated with 4 years' run at 
BEPCII. 



TABLE II: The sensitivity for Rm measurements at BES-III 
with different decay modes with 4 years' run at BESPCII 



D"W Mixing 





Reaction 




Events 


Sensitivity 








RS(xlO'') 


Rm {xlQ-^) 


i/'(3770) 


-> (A"7r+)(A'" 




10.4 


1.0 


V'(3770) - 


{K-e+v){K~ 


-e+v) 


8.9 




V'(3770) - 


{K-e+u){K- 




8.1 


3.7 


^/)(3770) - 


^ {K-^,+ v){K- 




7.3 





can obtain Rm = (1.18 ±0.6) x 10"". With 20fb-i data 
at BES-III, about 12 events for the precess D°D° 
{K^n^){K'^TT^) can be produced. One can observe 3.0 
events after considering the selection efficiency at BE- 
SIII, which could be about 25% for the four charged 
particles. The background contamination due to double 
particle misidentification is about 0.6 event with 20/5^^ 
data at BES-III j20i. Table |TTI] fists the expected mixing 
signal for N^.g = N{K^TT'^){K'^7r'^), background Nbkg , 
and the Poisson probability P{n), where n is the possible 
number of observed events in experiment. In Table IIIIl 
we assume the R]\i ~ 1.18 x 10^^, the expected number 
of mixing signal events are estimated with lOfb^^ and 
20fb^"'^, respectively. 



TABLE III: The expected mixing signal for Nsig — 
N{K^n^){K^TT^), background Ntkg , and the Poisson prob- 
ability P{n) in 10 fb"^ and 20 fb"^ at BES-III at V(3770) 
peak, respectively. Here, we take the mixing rate R]\f — 
1.18 X 10"". 





10 fb~^ (i/'(3770)) 
36 million D°D° 


20 fb-i (^/'(3770)) 
72 million D°D° 


Ns^g 


1.5 


3.0 




0.3 


0.6 


P{n = 0) 


15.7% 


2.5% 


P{n = 1) 


29.1% 


9.1% 


P{n = 2) 


26.9% 


16.9% 


P{n = 3) 


16.6% 


20.9% 


P{n = 4) 


7.7% 


19.3% 


P{n = 5) 


2.8% 


14.3% 


P{n = 6) 


0.9% 


8.8% 


P{n = 7) 


0.2% 


4.7% 


P{n = 8) 


0.1% 


2.2% 


P{n = 9) 


0.01% 


0.9% 



In the limit of CP conservation, by combing the mea- 
surements of X in KsTTTT and ycp from Belle, one 



In conclusion, we discuss the constraints on the strong 
phase difference in Ktt decay according to the 

most recent measurements of y' , ycp and x from B fac- 
tories. We estimate the sensitivity of the measurement of 
mixing parameter y at -0(3770) peak in BES-III experi- 
ment. With 20 fb~^ data, the uncertainty A(?/) could be 
0.003. Thus, assuming y at a percent level, we can make 
a measurement of y at a significance of 4.3(7 deviation 
from zero. The sensitivity of the strong phase differ- 
ence at BES-III are obtained by using data near the DD 
threshold with CP tag technique at BES-III experiment. 
Finally, we estimated the sensitivity of the measurements 
of the mixing rate Rm, and find that BES-III experiment 
may not be able to make a significant measurement of 
Rm with current luminosity by using coherent DD state 
at 0^(3770) peak. 
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